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Abstract 

We consider the differential conductance of a periodically driven system connected to infinite electrodes. We focus on 
the situation where the dissipation occurs predominantly in these electrodes. Using analytical arguments and a detailed 
numerical study we relate the differential conductances of such a system in two and three terminal geometries to the 
spectrum of quasi-energies of the Floquet operator. Moreover these differential conductances are found to provide an 
accurate probe of the existence of gaps in this quasi-energy spectrum, being quantized when topological edge states occur 
within these gaps. Our analysis opens the perspective to describe the intermediate time dynamics of driven mesoscopic 
conductors as topological Floquet filters. 
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1. Introduction 

Recently, the possibility to induce an out-of-equilibrium 
topological state of matter through irradiation or a peri¬ 
odic driving has stimulated numerous works. While ini¬ 
tially the external driving perturbation was used to trigger 
a phase transition between states of conventional topolog¬ 
ical order [1, 2, 3], fascinating topological properties spe¬ 
cific to driven out-of-equilibrium states were soon identi¬ 
fied [4, 5, 6]. While several proposals to realize and probe 
these topological states in various artificial systems have 
turned out to be successful [7, 8, 9, 10, 11, 12, 13, 14]. 
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Their realization in condensed matter have proved to be 
challenging [15, 16]. 

There is a strong analogy between equilibrium topolog¬ 
ical insulators and topological driven states. Both require 
the existence of a gap in the spectrum characterizing their 
single particle states: topological insulators are band in¬ 
sulators with a gap in the energy spectrum of the single 
particle Hamiltonian while topological driven states have 
a gap in the spectrum of the Floquet operator. In both 
cases, a nontrivial topology manifests itself through the 
appearance within this gap of robust states located at the 
edge of the system. However, while in an insulator the 
gap separates empty states from occupied states, the ther¬ 
modynamics of gapped periodically driven states is much 
less understood. Recent studies have stressed the differ¬ 
ences between the nature of the states reached at long 
time in such periodically driven systems and the equilib¬ 
rium ground states of insulators [17, 18, 19, 20]. 

Here we follow a different route: we focus on the re¬ 
lation between the DC transport of a periodically driven 
system and its quasi-energy Floquet spectrum in a regime 
where the times of flight of electrons through the sys¬ 
tem are shorter than the characteristic inelastic scatter¬ 
ing times, which can be the case in mesoscopic systems. 
This provides a way to avoid the issue of long time dy¬ 
namics of driven systems, which was raised in recent stud¬ 
ies [19, 20, 21]. Technically, this requires that the dom¬ 
inant perturbation of the unitary evolution of the driven 
system is the presence of the electrodes: dissipation should 
occur in the leads. From this point of view, the driven sys¬ 
tem behaves as a topological Floquet filter instead of an 
out-of-equilibrium steady-state analogue of an equilibrium 
insulating state. 

DC transport is known to be an ideal probe of the 
existence of edge states in topological equilibrium phases 
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realized in condensed matter, particularly in two dimen¬ 
sions. In a seminal paper [22] Markus Biittiker demon¬ 
strated how the non-local conductances in a Hall bar fully 
characterize the nature of the quantum Hall effect and the 
associated chiral edge states. This approach was recently 
extended to study the quantum spin Hall effect occuring in 
HgTe/CdTe quantum wells [23, 24]. In this time-reversal 
invariant topological phase, the existence of a Kramers 
pair of counter-propagating edge states leads to a series 
of non-local conductances whose experimental observation 
clearly identified this new phase. For topological driven 
systems, the situation is more confusing: building on ear¬ 
lier works on the transport through a topological period¬ 
ically driven state [25, 3], recent studies have focused on 
the transport through a one-dimensional topological su¬ 
perconducting state [26], the effect on transport of the 
competition between heating by the drive and the coupling 
to the leads [21] or the quantization of conductances of a 
topological phase in multi-terminal geometry [27]. It was 
also proposed to probe quasienergy spectra (and topologi¬ 
cal edge states) through magnetization measurements [28] 
and tunneling spectroscopy [29]. However, the relation 
between transport and the existence of topological edge 
states in periodically driven states remains unclear, and 
a summation procedure over different energies in the lead 
was proposed to recover a quantized conductance [26, 30]. 
The purpose of our paper is to reconsider the relation 
between the (non-local) differential conductances of pe¬ 
riodically driven systems and their Floquet quasi-energy 
spectrum, allowing for a direct relation between these dif¬ 
ferential conductances and the topological indices associ¬ 
ated with the spectral gaps. In particular we will establish 
a protocol in a multi-terminal geometry allowing for this 
identification. In this point of view, a topological periodi¬ 
cally driven system is viewed as a topological Floquet filter 
with selective edge transport occurring for specific voltage 
biases between a lead and the system. These voltage bi¬ 
ases lead to a stationnary DC current by counterbalancing 
the time dependence of Floquet states. 


constant in energy: ^{E) E. The dynamics of the sys¬ 
tem is described by the evolution operator U(t^t') which 
obeys the equation 

t') = - s) iJU t'). (1) 

Of great importance is the Floquet operator which is the 
evolution operator after one period t/(T, 0). When diag- 
onalizable, it can be decomposed on the left eigenstates 
{(j)cx\ and the right eigenstates \(j)a) of U{T^0) 

c/(r,o) 10a) = A„ !</>„), 

(4|C>(T,0) = Aa (0a I, 


that constitute a bi-orthonormal basis of the Hilbert space 

(0a I 0/3) =Sal3 ; Yi I = 

a 


The eigenvalues Aq, in Eq. (2) are called the Floquet mul- 
tiplicators and read 


Ack = exp 



- i7a) T 


(4) 


The coefficient is called the quasienergy and ja is its 
damping rate whose inverse gives the life-time of the eigen¬ 
state. Note that the quasienergy being a phase, it is de¬ 
fined modulo the driving frequency uj = 27r jT. Any state 
at arbitrary time t can then be constructed from the eigen¬ 
states of the Floquet operator. It is particularly useful to 
define the left and right Floquet states 


C/(i,O)|0a), 


which are periodic in time, \ua{t)) = \ua{t FT)) (same for 
(fta(^)l), so that they can be expanded in Fourier series 

(6) 


2. Prom Floquet theory to scattering theory 

2.1. Floquet theory for open systems 

We consider a periodically driven quantum system con¬ 
nected to Meads equilibrium electrodes through good con¬ 
tacts with large transmissions. The system is described by 
a Hamiltonian ^"^"(t) - E where H^^^t + T) = H^^^t) 
with T the period of the drive, and E is a self-energy ac¬ 
counting for the coupling between the system and its en¬ 
vironment (e.g. the leads). We assume in the following 
that this self-energy is dominated by the exchange with 
the electrons in the leads. When all characteristic times 
of the leads are small with respect to the characteristic 
times of the system, we can use the so-called wide band 
approximation [31] where the self-energy is assumed to be 


where the harmonics read 

\F) = ^ f 

^ (7) 

From Eqs. (3) and (5) the evolution operator can be ex¬ 
panded on the Eloquet states as 

Uit,t') = \u„{t)) {Ua{t')\. (8) 

a 

This expression can finally be decomposed on the harmon¬ 
ics of the Eloquet states by using Eq.(6) 

\u(p){uF\. (9) 

a 
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In practice, the spectrum of Floquet operator of the 
semi-infinite system can be obtained numerically either by 
direct a computation of U (T, 0) (e.g. as a discretized in 
time version of the infinite product) or through its repre¬ 
sentation in Sambe space [32]. 


which satisfy the rule Yhi (E) = 0 for any £' (the av¬ 
erage current leaving the system is independent on the 
chemical potentials in the leads). This formula is analo¬ 
gous to the Landauer-Biittiker formula for the differential 
conductance of multiterminal equilibrium systems. 


2.2. Differential conductanee 

Based on a standard formalism, we can calculate an¬ 
alytically the differential conductance of the periodically 
driven system in a multi-terminal geometry and relate it 
to the quasienergy spectrum of the system. We follow the 
standard Floquet scattering formalism [33, 31, 34, 35, 26, 
27] to describe the transport properties of this multiter¬ 
minal setup in a phase coherent regime. We consider the 
rolling average over a period T (all time-averages in the fol¬ 
lowing are also rolling averages over one driving perdiod) 
of the current entering each lead labelled by the index t. 

= {Je{t')). (10) 

where is the expectation value of the current en¬ 

tering lead £ at time t'. This average current satisfies a 
relation [33, 31, 34]: 


2.3. Generalized Fisher-Lee relations 

The average transmission coefficients T££f{E) can be 
related to the Floquet-Green functions of the system, in 
a way analogous to the case of undriven conductors [36]. 
The retarded Green functions in mixed energy-time repre¬ 
sentations are defined by 

Q{t,E) = ^ J dr g{t, t-r) (15) 

where Q{tD') satisfies 

git, t') = 5it - t'). (16) 

It is related to the evolution operator defined in Eq.(l) by 
git,t') = fQit-t')Uit,t'). (17) 


liit) = ^ dEYi [Tee>it,E)fe,iE)-Te,eit,E)feiE)], 

( 11 ) 

where f£{E) is the Fermi-Dirac distribution of the lead 
£ assumed to be at equilibrium at the chemical potential 
ja£. The T££>{E) are the time-averaged transmission coeffi¬ 
cients between lead £ and £' which will be discussed below. 
We define the differential conductance G££'{t^E) as the 
sensitivity of the current entering the lead £ to variations 
of the chemical potential /i^/ of the lead^ £' : 

Gu'it,E) = -e^ . ( 12 ) 

Note that this definition is not symmetric in the vari¬ 
ous chemical potentials as opposed to other defini¬ 
tions used in the literature. In the long time stationary 
regime on which we focus, t ^ oo, the average conduc¬ 
tances G££'{t^E) are expected to reach a value indepen¬ 
dent on the choice of origin of time t and associated initial 
conditions [31], which we denote by G^,{E). 

We obtain from Eq. (11) the zero temperature time- 
averaged differential conductances 


Erom the decomposition Eq.(9) of U{tD') over the har¬ 
monics of the Eloquet states, we express the Floquet-Green 
funetions 


= ^ r dte^p^^gUE). 

^ Jo 


as 


gip)iE) = Y 




E — [Sa F rhuj — ih^a] ' 
The transmission coefficients are expressed as 


Te>eiE) = YTiFE), 


(18) 

(19) 


( 20 ) 


where each T^fj{E) is the transmission coefficient for an 
electron injected in lead £ at the energy E = fi£ and leaving 
the system in lead F at the energy E -\- phuj^ i.e. after 
having exchanged p quanta with the driving perturbation. 
They read (see Appendix A and [31, 34, 35, 37]): 


T^f{E) = Tr g^^^\E)t£>{EFphuj)g^^\E)t£{E) (21) 


Gf,, (E) = jTu> (E) for i ^ I' (13) 

GT,iE) = -'£-YniiE), (14) 

^See [35] for a recent discussion of chemical versus electrical po¬ 
tential drops at the interface between the system and an electrode. 
Such subtleties are quietly neglected in the following. 


where f^(E) is the coupling operator at energy E between 
the system and the electrode £. 

Note that when the rates ja are sufficiently small for 
the quasi-energies occurring in Eq. (19) to be “well- 
defined”, the equations (19,20,21) imply that the transmis¬ 
sion coefficient T£>£{E)) and thus the differential conduc¬ 
tance G£>£{E) vanishes whenever the energy ii£ = E does 
not correspond to a quasi-energy e^ up to a multiple of huu. 
This is nothing but the conservation of energy of incident 
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Figure 1: Lattice of couplings for the half-BHZ model with two 
orbitals per site of the square lattice. Blue and red links correspond 
to the and nearest neighbor coupling matrices between s 

and p orbitals while the green links represent the second-nearest 
neighbor coupling matrix. 


state which holds modulo huj in a periodically driven sys¬ 
tem. This property allows one to probe the existence of 
gaps in the quasi-energy spectrum of the Floquet opera¬ 
tor. Note however that to access the whole quasi-energy 
spectrum, the lead jui has to be strongly biased (with a 
bias of order huj) with respect to the scattering region at 
/isys, a situation far from the standard procedure to probe 
equilibrium phases, but required to probe an inherently 
out-of-equilibrium Floquet state. Moreover, in the pres¬ 
ence of a ballistic mode connecting two leads £ and £' such 
as the edge mode of a topological Floquet state, we expect 
the corresponding Tifi{E) to be quantized provided both in¬ 
terfaces are sufficiently transparent. Note however that as 
the states leaving the systems have component of energies 
E phw, this implies a perfect transparency over a broad 
range of energies. These two properties allow for a poten¬ 
tial probe of topological Floquet states through non-local 
transport. In the next section, we study this application 
by a numerical implementation of transport through an 
AC driven system. 



IF 
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Figure 3: Geometry of the semi-infinite system used to compute 
quasi-energy spectra of figure 2. In those spectra, the color of each 
point represents the average y position of the corresponding Floquet 
mode, according to the color scale shown on this figure. 


The tight-binding Hamiltonian with nearest and next-to- 
nearest neighbors hoppings (see Fig. 1) can therefore be 
written as a two-by-two matrix on the (s,p) basis as 


^BHZ = 


E 

x,y 


[{M - J - 4B) - jiao] |a;, y) (x, yj 


+ tP \x,y) (a; + l,y| \x,y) {x,y+l\ 

+ {\x,y) {x + l,y + l\ + \x,y) {x + l,y - l\) 


( 22 ) 


-f h.c. 


with hopping matrices 

and = (23) 

where (7x,y,z are the Pauli matrices and ao the identity 
matrix. The parameters of this Hamiltonian are chosen 
as M = —1.0, J = 1.5, A = 4.0, B = 1.5 : they correspond 
to an equilibrium phase which is a trivial insulator. A 
topological insulating phase can be reached by varying e.g. 
M such that 0 < M/B < 4 or 4 < M/B < 8. This 
trivial equilibrium phases is submitted to a periodic on¬ 
site perturbation [5] 

Ai7(t) = E [sm{ujt)cFx + cos{ujt)ay] |x, y) {x,y \ . (24) 
x,y 


3. Numerical study 

In the following, we provide numerical evidence of the 
correspondence between the differential conductance and 
the quasi-energy spectrum in a periodically driven system. 
In particular, we show that a single out-of-equilibrium 
topological edge state corresponds to a quantized differen¬ 
tial conductance. The chirality of these topological edge 
states is probed in a multi-terminal setup. 

3.1. Model and method 

Following [5], we use the restriction to spins up of 
the Bernevig-Hugues-Zhang model of quantum spin Hall 
effect in HgTe/CdTe quantum wells [38] (referred to as 
half-BHZ model). As the Haldane model [39] it realizes 
an anomalous quantum Hall equilibrium phase, but on a 
square lattice with two orbitals per site, denoted s and p. 


Note that this perturbation does not correspond to the 
variation of a parameter of the initial Hamiltonian. Through¬ 
out our study, we have used a driving frequency cj = 20. 
Depending on the strength F, this perturbation can drive 
the system either towards a topologically nontrivial out- 
of-equilibrium state or in a topologically trivial out-of¬ 
equilibrium state. We have chosen two values E = 2 
and F = 8 of the driving amplitude so that gaps open at 
e = 0, Ehuj/2 of the quasi-energy spectrum. For each value 
of F and each gap, we have computed numerically the bulk 
topological invariant IFe[t/] associated to the quasi-energy 
gap around e [5] to ensure that F = 2 correspond to a 
trivial gapped Floquet state, while F = 8 corresponds to 
a topological gapped Floquet state with a non-trivial gap 
at e = 0. Alternatively, in the two cases we have computed 
the quasi-energy spectrum for the driven model in an in¬ 
finitely long ribbon of width IF = 60 sites (see Fig. 3) by 
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Figure 2: Quasi-energy spectrum of the driven half-BHZ model in an infinitely long ribbon and corresponding asymptotic two terminal 
differential conductance G^-^(E) for finite system connected to infinite leads. The results are shown for parameters with only trivial gaps 
with driving amplitude F = 2 (Left) and a topological gap at e = 0 or F = 8 (Right). The quasi-energy spectra are obtained for strip of 
width VF = 60 sites by diagonalization in Sambe space with 5 sidebands. The quasi-energies e are shown in units of the driving frequency 
huj as a function of the quasi-momentum in the direction of the ribbon kx- For the topological case (F = 8.0, Right figure), states appear 
inside the gap e = 0, located on either side of the ribbon. The location of each Floquet mode was numerically computed and converted to a 
color using the color scale of figure 3. Delocalized states (and states located in the bulk) are in green, whereas states located on an edge are 
in red or blue, depending on which edge. This signals the appearance of a topologically non-trivial gap around e = 0. On the other hand, the 
gaps in e = remains topologically trivial in both cases. The asymptotic long time limit of the corresponding two-terminal differential 

conductance is displayed as a function of the chemical potential of the left lead (black dots and lines) for a finite sample of size VF x L = 60 x 30 
sites. This differential conductance quantized, with G = \e^ jh in the topological quasi-energy gap (Right) and vanishes in the trivial 

gap (Left). Hence it is an accurate probe of the existence of topological gap in the quasi-energy spectrum of open conductors connected to 
electrodes. The colored arrows indicate the chemical potential of the incoming lead used in the computations of Fig. 4. 


diagonalization in Sambe space with 5 sidebands. The re¬ 
sulting spectra are shown in Fig. 2: they show as expected 
that states located at each edge on the ribbon appear in¬ 
side the e = 0 gap for the topological case (F = 8) as 
opposed to the trivial case (F = 2). 

Finally, to study transport through the system, leads 
are attached to a finite size system. These leads are mod¬ 
eled by a simple tight-binding Hamiltonian on a square lat¬ 
tice with nearest neignbors hoppings with amplitude Jo = 
8. An onsite potential is added to the lead Hamiltonian 
to reduce mismatches between the incoming and outgoing 
states of the leads and the scattering states of the central 
region. 

The numerical calculations are performed using the nu¬ 
merical method described in [35]. Although the technique 
is based on wavefunctions, it is mathematically equivalent 
to the Green function approach used in this article (see 
Eq. (49) of [35] for the connection with the transmission 
coefficient as well as section 5.4 for the link with Floquet 
theory). Our implementation is based on the Kwant pack¬ 
age [40]. 

3.2. Probing the quasi-energy bands and gaps through two- 
terminal differential conductance 

The differential conductance GRL(t) is computed in a 
two-terminal setup through a sample of width W = 60 and 
length F = 30 sites. Throughout this study, the chemical 
potential of the undriven system was chosen as //gys = 0- 
To get rid of oscillations faster than the driving frequency. 


which are irrelevant for our purpose, we perform a slid¬ 
ing average over one driving period T. After a transient 
regime, the (averaged) differential conductance G(t) con¬ 
verges to a finite value (see Fig. 4). This transient regime 
can be understood as the time of flight of the state in¬ 
jected at the left lead to the right lead after the driving 
perturbation has been turned on and the Floquet states 
developed inside the system. When the chemical potential 
of the incoming lead lies in a topological quasi-energy gap 
of the scatterng region, transport occurs through a chiral 
state localized near the edge of the sample. We can easily 
evaluate the travel length L. The expected group velocity 
is extracted from the slope of the quasi-energy dispersion 
relation from Fig. 2 through 


and we obtain = 0.15(1) acj where a is the lattice spac¬ 
ing and uj the driving frequency. This correlates perfectly 
with the time At of the first increase from zero of the con¬ 
ductance from the switching on of the driving field. From 
the curve on Fig. 4 we find find F/At = 0.14(1) auj. 

After this transient regime, the differential conductance 
GRL(t) reaches a long time stationary limit G^l, as shown 
in the inset of Fig. 4. As expected, this asymptotic differ¬ 
ential conductance is sensitive to the quasi-energy spec¬ 
trum of the driven system: when in a spectral band, it 
reaches high values but vanishes when the chemical poten¬ 
tial /iL — Fsys lies in a trivial spectral gap of the Floquet 
operator, as shown on the left plot of Fig. 4. Moreover, 
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Figure 4: Two-terminal differential conductance as a function of time, in the transient regime, for two chemical potentials, in a topologically 
trivial case. Green curve: incoming chemical potential in a quasi-energy band ; blue and red curves: incoming chemical potential in the 
bulk quasi-energy gap, in a trivial (blue) and a topological (red) case. Inset: longer simulations were carried out to ensure convergence. The 
asymptotic value is denoted by 


when this spectral gap is topological, the associated pres¬ 
ence of chiral states at the edge of the system shown in 
Fig. 2 (Right) leads to a quantized two terminal conduc¬ 
tance as shown in the right plot of Fig. 4. To further study 
the correlation between the differential conductance as a 
function of the chemical potential /ip — Msys = E and the 
spectral gap of the system we have studied the behavior 
of the long time limit of this conductance as a function 
of E. The corresponding results as well the spectra of the 
isolated driven models are plotted in Fig. 2 for both the 
trivial and a topological gapped Floquet states. We find a 
strong correlation between the vanishing of both quantities 
for the topologically trivial case: the differential conduc¬ 
tance vanishes only inside a trivial gap, except at the edge 
of the gap where finite size effects occurs due to imperfect 
transparencies of the contact with the leads, as shown on 
a magnified view around the gap e = 0 in Fig. 5. This 
demonstrates both that the quasi-energy spectrum of the 
finite system connected to infinite leads is sufficiently close 
to the spectrum of the isolated infinite strip, and that the 
differential conductance is an accurate probe of this spec¬ 
trum for the open system. Moreover, in the topological 
case the asymptotic differential conductance remains con¬ 
stant equal to the number n = 1 of edge states (in units 
of e^/h) inside the topological gap as shown in Fig. 5. 
The small deviations from /h visible in fig¬ 

ure 5 are attributed to the finite dispersion relation of the 
leads, which does not completely satisfy the wide band ap¬ 
proximation and leads to imperfect transparencies of the 
contacts. 

3.3. Multiterminal geometry 

A crucial characteristics of the edge states occurring 
in the topological gap of the spectrum shown in Fig. 2 
(Right) is their chiral nature. In the equilibrium case this 



Figure 5: Zoom of Fig. 2 : asymptotic two-terminal differential con¬ 
ductance units of e^//i with respect to the chemical potential 

of the left lead (black dots and lines) around the quasi-energy gap 
e = 0 for a trivial (red) and topological gap (red). The reference 
energy /j,i^ = 0 corresponds to an unbiased lead whose chemical po¬ 
tential identifies with the one of the driven system /igys = 0. The 
differential conductance vanishes in the quasi-energy gaps for a triv¬ 
ial spectral gap, but reaches an approximately constant and quan¬ 
tized value G = e^/h within the topological gap with one chiral edge 
mode. Deviations from Gg^ = le^/h are attributed to the finite 
dispersion relation of the leads, which does not completely satisfy 
the wide band approximation and leads to imperfect transparencies 
of the contacts. 

property together with their ballistic propagation leads to 
quantized conductance in a Hall bar geometry [39]. To test 
the chirality of the topological edge states, we have com¬ 
puted differential conductances in a three-terminal geom¬ 
etry shown in Fig. 6, where the width of the contact with 
the electrodes is IF = 30 sites, the total length (between L 
and R contacts) is 50 sites, corresponds to all three arms 
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having a length of 10 sites. In this geometry, we monitor 
the two differential conductances 


Gr,t{E) 




<ELr 


^^T=|~lsys+E 


dlL 

dfiT 


flT —Msys^r E 


(26) 


where are the current in the i?, L contacts averaged 
over one period T of drive (see eq. 10). The chemical 
potential of the system is still set to /igys = 0 (as if imposed 
e.g. by a backgate). 
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Figure 6: Three-terminal differential conductance Glt (violet) and 
(orange) as a function of time for a chemical potential of the 
top lead in a topological bulk gap. After an initial transient regime, 
asymptotic values for the nonlocal conductances are reached with 
^LT ~ 0.0002(1) e^/h and = 0.94(1) e^/h. The simultane¬ 
ous vanishing of and quantized to the two-terminal value 

demonstrates the chiral nature of topological edge states. This is 
further illustrated in the inset by the corresponding probability den¬ 
sity |'0(t, a:)p in the stationnary regime (here at t = 1000 T) which 
is entirely localized around the T ^ R edge, demonstrating further¬ 
more the high transparency of the R contact for this edge mode. 
Simulations are done with a width of the contacts IT = 30 sites and 
a total length of 50 sites. 

We consider the case of a topological gap round e = 0 
(see Fig. 2 (Right)). First we set the chemical potential of 
the top lead inside this gap {rt = 0). The time evolution 
of the differential conductances are shown in Fig. 6. After 
a transient regime, the differential conductances converge 
to asymptotic values = 0.0002(1) e^/h and = 

0.94(1) e^/h. The value of G^ is in agreement with the 
two terminal results, while the vanishing of G^ is in per¬ 
fect agreement with the chiral nature of the edge state 
moving clockwise for the chosen parameters. This can be 
contrasted with the case of a chemical potential /it = 0 
set inside a bulk Floquet band, displayed in Fig. 7. In this 
case, after a longer transient regime due to slower group 
velocities, both conductances converge towards large finite 
values, confirming the absence of chirality for these bulk 
states. 

Finally, to illustrate the spatial structure of these Flo¬ 
quet states we show in insets of Fig. 6 and 7 a color map of 



0 500 1000 1500 2000 

tjT 


Figure 7: Three-terminal differential conductance Glt (violet) and 
Grt (orange) as a function of time for a chemical potential of the top 
lead in a quasi-energy band. A longer transient regime is observed, in 
agreement with smaller group velocities than for the edge states, and 
large asymptotic values for both Glt and Grt are reached, confirm¬ 
ing the non-chiral nature of the corresponding Floquet states. The 
inset shows color map representing the probability density 
for at t/T = 1800, demonstrating the delocalized nature of the Flo¬ 
quet states. 


the probability density of the states reached at large times 
for both the topological gap and inside the bulk bands. 
The chiral nature of the states injected in the top lead is 
clearly apparent and confirms the transport results. 


4. Discussion 

We have studied the differential conductances Gi>i{E) 
defined in Eq. (12) for both a two-terminal and a three- 
terminal geometry. As for equilibrium phases, we have 
found that these conductances probe the topological na¬ 
ture of a gapped Floquet state: it vanishes whenever the 
chemical potential ri lies in a topologically trivial gap 
while we found it to be quantized for chemical potential 
inside a topological gap of quasi-energies. Moreover, the 
chiral nature of the edge states reflected on the value of the 
multi-terminal differential conductances within this gap. 
These results validate the relation between the DC trans¬ 
port through a periodically driven system and the Flo¬ 
quet description of the dynamics of the isolated system. 
This requires that the dissipation occurs mostly inside the 
leads, i.e. that the system is small enough for the travel 
times through it to be small compared with times scales 
of other source of dissipation (phonons, photons, etc). In 
that situation, the differential conductance of the system 
depends on the intermediate time evolution of the driven 
system and can be accurately described starting from the 
unitary dynamics of the driven system and its associated 
topological properties [5, 6]. This should be the case as 
long as phase coherence is preserved on the length scale 
of the sample, even in the presence of dissipation or weak 
interactions. This result paves the way towards a direct 
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engineering of a topological filter through the implemen¬ 
tation of Floquet states. 

Our results on the quantization of differential conduc¬ 
tances inside a topological gap contrast with some of the 
recent results on similar systems [26, 27, 30], where a sum¬ 
mation procedure over multiple chemical potentials was 
required to recover a quantized conductance. Note that 
these previous results used a different measuring scheme 
of differential conductance, symmetric in chemical poten¬ 
tials of electrodes, which may filter out some of the energies 
of the outgoing electronic states, and thus require a sum¬ 
mation protocol. In our case, a strong and asymmetric 
voltage bias is used to compensate the time dependence 
(through quasi-energies) of the Floquet states. Moreover, 
we have used in our numerical study a model with a single 
edge mode located around k 0 diS opposed to double¬ 
mode topological gaps used previously. More complex edge 
structures may well lead to lower transparency of the inter¬ 
faces for energies in the gap, and thus non-quantized con¬ 
ductances. The study of the transparency of such interface 
between a strongly driven conductor and a DC electrode is 
a deep and topical subject whose discussion goes beyond 
the scope of the present study. 
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Appendix A. Transmission coefficients 

In the following, we extend the approach of [31] to 
a quasi-unidimensional system in the geometry of Fig. 3. 
The purpose is to describe the scattering through a two di¬ 
mensional driven system, viewed as a Chern Floquet filter. 
We follow closely the derivation of [31] and only highlight 
the differences due to the geometry considered. 

We consider a system composed of a central region IZ 
(the scattering region), which is submitted to a periodic 
excitation. This region is described by a time-periodic 
(with period T = 27r/co’) tight-binding Hamiltonian 

Hsys = Hnn' (^- 1 ) 

nn'£lZ 

with + T) = where is the annihilation 

operator of an electron in an localized state n of the tight- 
binding model, the index n representing the position on 


the Bravais lattice as well as internal degrees of freedom 
(sublattice, orbital, spin, etc.). This central region is con¬ 
nected to Meads leads described by the Hamiltonian 

-Nleads 

-f^leads ~ EE 

r=l qoc 

where a labels transverse modes of the semi-infinite lead 
and q is the longitudinal momentum in this lead. We as¬ 
sume that the wavefunctions in the leads read 

'lpag{x,y) = (A. 3 ) 

and in particular that transverse modes |x) do not depend 
on q and constitute an orthonormal basis of the trans¬ 
verse Hilbert space,. (Notice that this is not always the 
case, especially when a magnetic field is present [40]). The 
annihilation operator of a state localized at transverse po¬ 
sition y = yi{n) in the interface with the lead ^ is written 
as 

^i,y=yi{n) ^ ^ ^£qaXia (^) • (^•^) 

qa 

The internal degrees of freedom in the leads can be taken 
into account, if needed, by considering more leads. 

The contacts between the central region and the leads 
are described by the Hamiltonian 

-Nleads 

^contacts = ^ ^ ^ ^ (^•^) 

£=l neMi q 

where Mi describes the set of sites of the central region at 
the interface with lead £. In terms of the transverse modes 
creation/annihilation operators, this Hamiltonian reads 

Weads 

ffcontacts= E E WMcn + h-C.. (A. 6 ) 

£=l nEA4i qa 

Using the approach of [31] which amounts to describe 
the correlations of the driven system in terms of the equi¬ 
librium noise in the leads, we obtain the equation (11) of 
main text with transmission coefficients 

yp{E) = 

m',n'GX£/ 


that can be written as a trace on the interfaces (Eq. (21) 
of main text). 


T^l){E) = Tr \[g^P'>{E)£T\E + phuj)g^P'>{E)T^'{E) 


where 

TeiE) = 27 r^ |V,|" \xa){Xa\S{E - E^,^). (A. 7 ) 

qa 

This expression is in agreement with previous approaches 
based on different formalisms [34, 35, 37]. 
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